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Abstract. 

We have re-examined an ancient VLBI survey of ultra-compact radio sources at 
2.29 GHz, which gave fringe amplitudes for 917 such objects with total flux density 
> 0.5 Jy. A number of cosmological investigations based upon this survey have been 
published in recent years. We have updated the sample with respect to both redshift 
and radio information, and now have full data for 613 objects, significantly larger 
than the number (337) used in earlier investigations. The corresponding angular- 
sizc/redshift diagram gives Q m = 0.25 + 0.04/- 0.03, n A = 0.97 + 0.09/- 0.13 and 
K = 0.22 + 0.07/— 0.10 (68% confidence limits). In combination with supernova data, 
and a simple-minded approach to CMB data based upon the angular size of the acoustic 
horizon, our best figures are Q m = 0.304 + 0.024/- 0.023, fl A = 0.693 + 0.034/- 0.035 
and K — —0.003 + 0.021/— 0.019. We have examined a simple model of vacuum energy, 
based upon a scalar potential V(<p) — loq<\> 2 /2, to test the possibility that the vacuum 
is dynamical; whereas the data favour a value ujq = 0, they are compatible with values 
in excess of the Hubble rate. 



1. Introduction 



The angular-size/redshift relationship is in principle a simple cosmological test, which 
has yet to play a noted part in the development of observational cosmology. Endeavours 
in this field have generally concentrated upon classical double radio sources as putative 
standard measuring rods, angular size being defined by core-lobe or lobe-lobe separation. 
Early results were problematical, showing angular sizes which diminish more rapidly 
with increasing redshift z than would be allowed by any plausible cosmological model. 
This behaviour is believed to be an evolutionary effect, brought about by interaction 
with an evolving extra-galactic medium (Legg 1970; Miley 1971; Barthel and Miley 1988; 
Singal 1988), or a selection effect, due to an inverse correlation between linear size and 
radio power (Jackson 1973; Richter 1973; Masson 1980; Nilsson et al 1993). There have 
been several attempts to disentangle the latter effects from proper cosmological ones 
(Daly 1994; Buchalter et al 1998). Daly (1994) introduces a model of lobe propagation 
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which allows intrinsic size to be estimated from other source parameters; Buchalter et 
al (1998) select Fanaroff- Riley Type II sources (Fanaroff and Riley 1974) with z > 0.3 
as the basis of their work, in part because their morphology allows an unambiguous 
definition of angular size, and also allow for a linear-size/redshift correlation of the sort 
discussed above; such measures bring the angular-size/redshift diagram for extended 
sources into concordance with acceptable Friedmann cosmological models, but not with 
sufficient precision to allow definitive statements about cosmological parameters. 

More promising candidates in this context are ultra-compact radio sources, with 
milliarcsecond angular sizes measured by very-long-baseline interferometry (VLBI) 
(Kellermann 1993; Gurvits 1994). Gurvits work was based upon a large VLBI 2.29 GHz 
survey undertaken by Preston et al (1985) (hereafter referred to as P85). The latter was 
designed to establish a comprehensive full-sky list of ultra-compact radio sources; 1398 
candidates were selected from existing surveys, on the basis of total flux density at 2.7 
GHz. Selection criteria were complicated, but the candidate list is believed for example 
to be 97% complete to a limit of 1.0 Jy for sources with spectral index a > —0.5, 
which comprise about 80% of the list; further details are given in the original reference. 
These sources were then observed in a series of VLBI experiments at 2.29 GHz; compact 
structure was detected in 917 cases, and the corresponding correlated flux density S c 
(fringe amplitude) recorded, plus the total flux density S t if this was measured at the 
same time (531 out of the 917 cases). As angular size is defined by fringe visibility 
S c /S t (Thompson et al 1986), the number of potential candidates for an investigation 
based upon P85 is thus appears to be 531. In relation to these P85 give 269 redshifts; 
Gurvits (1994) added further redshifts from Veron-Cetty and Veron (1991), to give a 
sample of 337 objects in all. Gurvits gave plausible reasons for ignoring sources with 
z < 0.5, and using just the high-redshift data (258 objects) found marginal support 
for a low-density Friedman cosmological model, but considered only those with vacuum 
energy density parameter JIa = 0. Using exactly the same data set (kindly supplied by 
Dr. Gurvits), Jackson and Dodgson (1997, submitted 3/05/96) extended this work to 
the full density-parameter plane, and concluded that the best values are Q m = 0.2 and 
tt A = 0.8, if the Universe is spatially flat, later refined to fi m = 0.24 + 0.09/ -0.07 (95% 
confidence limits) (Jackson 2004). 

The latter work develops an astrophysical model, according to which the underlying 
population consists of compact symmetric objects (Wilkinson et al 1994), comprising 
central low-luminosity cores straddled by two mini-lobes. Ultra-compact objects are 
identified as cases in which the lobes are moving relativistically and are close to the line 
of sight, when Doppler boosting allows just that component which is moving towards the 
observer to be observed; the interferometric angular sizes upon which this work is based 
correspond to the said components. As z increases a larger Doppler factor is required; 
it turns out that latter approximately cancels the cosmological redshift, so that the 
observed component is seen in its rest frame. This model, coupled with the fact that 
the central engines which power these sources are reasonably standard objects (black 
holes with masses close to 1.5 x 10 10 M Q ), gives a plausible account of their behaviour as 
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standard measuring rods, and of why those with z < 0.5 should ignored. Those having 
z > 0.5 are the ultra-compact analogues (and perhaps the precursors) of Fanaroff-Riley 
Type II sources, and there are parallels between this work and that of Buchalter et al 
(1998). This is a slightly oversimplified precis of the account given in Jackson (2004), 
to which the reader is referred for further details. 

We have taken a fresh look at the P85 catalogue, and updated the latter with respect 
to both redshift and radio information. Details of the new data set are given in Section 
2. Section 2 also discusses the acoustic horizon and its linear size; coupled with accurate 
measures of the location of the first peak in angular spectrum of the cosmic microwave 
background (CMB) radiation (Balbi et al 2000; de Bernardis et al 2000; Hanany et 
al 2000, Hinshaw et al 2003), these allow the WMAP results to be interpreted in a 
simple and transparent fashion, which does not need the full machinery of CMBFAST 
(Seljak and Zaldarriaga 1996), that is if particular figures for Hubble's constant and 
baryon content are accepted. Cosmological results are presented in Section 3, for the 
new data alone and in combination with the most recent supernova observations (Riess 
et al 2004; Astier et al 2006), mainly as marginalized confidence regions in the f2 m -^A 
plane. Section 4 considers the evolution of vacuum energy in the light of the new data. 

Although we believe that the results to be presented here are of intrinsic interest, a 
large part of our motivation in undertaking this work is to enhance further the credibility 
of milliarcsecond radio-sources as standard measuring rods. In this context it should be 
noted that their first application to the fi m -f2A problem (Jackson and Dodgson 1997) 
was blind, in the sense that it gave answers close to the concordance model at a time 
when Type la supernovae were giving fi m ~ 1 (Perlmutter et al 1997). This situation did 
not change until the advent of Schmidt et al (1998), Riess et al (1998) and Perlmutter 
et al (1999). 

2. Data 

Of the 1398 sources in the P85 catalogue, 917 have definite measures of correlated 
flux density (the rest have upper limits), and 531 have both total and correlated 
flux densities. We have consulted the NASA/IPAC Extragalactic Database (NED), 
and find 957, 723 and 456 redshifts for the respective sub-samples (end of August 
2005). Thus we have redshifts and visibilities for 456 sources; formally 8 of these 
have visibility>l and were discarded; the remaining 448 comprise our gold standard 
sub-sample. There are 386 sources for which a correlated flux density is listed but no 
total flux density (531+386=917); we have attempted to update this sub-sample, by 
searching for convenient listings in the literature which are roughly contemporaneous 
with P85 and at a frequency not too far removed from 2.29 GHz. With respect to 
southern hemisphere sources, we find that the Parkes catalogue PKSCAT90 (Parkes 
Catalogue 1990) serves this purpose very well. The latter gives total flux densities at 
2.7 Ghz, which we have extrapolated to 2.29 GHz using the spectral indices given in 
P85 (in the few cases where a is not listed we have used a value of zero). This procedure 
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gives 256 new values of S t , for which we have 174 redshifts; formally 9 of these have 
visibility>l and were discarded; the remaining 165 sources comprise our silver standard 
sub-sample, which we expect to be of lower quality than the gold. We have checked 
the validity of this approach by considering the overlap: Figure [T] is a plot of the total 
flux density so extrapolated against that listed in P85, for all sources (346) for which 
we have both; the two figures are well correlated, with a correlation coefficient of 0.854. 
With respect to northern hemisphere sources, we have examined the NRAO/VLA Sky 
Survey (NVSS) (Condon et al 1998), which gives St at 1.4 GHz, but find that this data 
set (our bronze sub-sample) is less useful, being too far removed in frequency and epoch. 
The sample used here will comprise our gold plus silver sub-samples, amounting to 613 
objects out of a possible 917, with a range z = 0.0035 to z — 3.787. 

Preston et al (1985) selected their list of 1398 VLBI candidates mainly from an 
earlier version of the Parkes survey (PKSCAT85, see Parkes Catalogue 1990), and 
from the NRAO-Bonn survey (Kuhr et al 1981); it would have been useful for this 
investigation if the corresponding flux densities had been included in P85, but this was 
not the case. Almost all of those sources in Kuhr et al (1981) were re-measured at 2290 
GHz by Preston et al (1985). 

Regarding the acoustic horizon, there are several publications which give similar 
analytical formulae for its linear size /acH (Hu and Sugiyama 1995; Mukhanov 2004). 
Writing that given by Hu and Sugiyama (1995) in terms of quantities which are all 
evaluated at recombination of the photon-baryon plasma, we find 



iidI =ViI ln l aTTTVi — (1) 



where 



A= = [±*Ll^-\ (2) 

v (pmY 1 ( Pc M m h 2 \ ~* 

C = 2c {^^)^ = 2c(Si n Xr"*(l + z,r 3 ' 2 - (4) 

Here p c is the critical density corresponding to a particular value of Hubble's constant 
Hq = 100/i km sec -1 Mpc -1 ; subscripts b and m refer to baryons and baryons+cold 
dark matter, and 7 and R to CMB photons and photons+neutrinos, respectively; a 
subscript r refers to values at recombination, and zero to current values. Equation (0) 
is exact only in the case of flat cosmologies, but in practice flatness will hold to a high 
degree of accuracy for any realistic model, over the interval between the big bang and 
recombination. It is exact also only in the case of instantaneous recombination at some 
redshift z T , which is unrealistic; in reality recombination starts at a redshift of about 
1200 and is complete at about 800. For this reason equation (JTJ) underestimates the 
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size of the acoustic horizon, because it does not allow for those charged baryons which 
recombine in the period immediately preceding z r , and curtails acoustic propagation 
too early (Mukhanov 2004). We have allowed for this by devising a semi-empirical 
two-parameter fitting formula, which increases /acH by a factor a(Qb,h), and reduces 
Pr by a fixed factor (3\ the latter is effective because it allows more time for acoustic 
propagation to take place. Taking z r = 1088 (Hinshaw et al 2003) and three massless 
neutrinos, we find that the following values are appropriate § : 

a = 1.2047 + 1.008(O b - 0.04) + 0.1178(/i - 0.7) (5) 

(3 = 0.602. (6) 



Comparison with numerical results from CMBFAST show that with these modifications 
equation (JTJ) is accurate to better than 0.4% over the ranges 0.03 < Q\> < 0.05, 
0.65 < h < 0.75 and 0.2 < Q m < 0.5, and to better than 0.6% if the latter range 
is increased to 0.2 < Q m < 1.0. We convert ZacH into an angular scale using the correct 
angular-diameter distance cIa(z), which need not correspond to a spatially flat model. 
As an example, Figure 121 shows the multipole index Z pea k corresponding to the first 
Doppler peak in the CMB angular spectrum as a function of f2 m , according to equation 
in the case of spatial flatness, Qa — 1 — ^m- 

For Type la supernovae (SNe la) we have used a composite data set comprising 
objects from Riess et al (2004) (the gold sample) and Astier et al (2006) (The Supernova 
Legacy Survey, SNLS). Each of these samples includes a low- 2 subset (z < 0.125, 68 and 
44 objects respectively); the two low-z subsets are not identical, but have 21 objects 
in common. Our composite data set comprises the full gold sample from Riess et al 
(2004) (157 objects) plus 71 high-z objects from Astier et al (2006) (Table 9 of their 
paper); to make sure that no objects were counted twice, we do not include the SNLS 
low-;? subset in our composite sample. However, before assembling the latter there is a 
matching problem to be considered, discussed below. 

The distance indicator listed in both Riess et al (2004) and Astier et al (2006) 
is distance modulus fiB', to bring the two samples into conformity we find that it is 
necessary reduce the gold subset distance moduli by 0.14 magnitudes. This figure was 
determined by fitting the equation 

fi B = 5log 10 (z) + C (7) 

first to the low-z gold subset, and then to the low-z SNLS one; the best values are 
C = 43.37 and C = 43.23 respectively. We choose the low- z subsets for this purpose 
(rather than the full ones), because they overlap, and because the outcome does not 
depend upon a choice of cosmological parameters. 

The probable reasons for this offset are outlined below. Astier et al (2006) list a 
rest-frame maximum-luminosity magnitude m* B for each supernova, from which they 

§ Including z r dependence the equation is: a = 1.1991 + 7.0 x 10~ 4 (,z r — 1080) 

+ [1.020 - 1.5 x 10- 3 (z r - 1080)](Q b - 0.04) + [0.121 - 4.0 x 10~ 4 (z r - 1080)] (/i - 0.7) 
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derive a corresponding magnitude m# corrected for brighter-slower and brighter-bluer 
correlations, which reduce the dispersion in magnitude by a factor of about 2 (see 
for example Guy et al 2005 and the references given there). The magnitudes so 
corrected were used in a Hubble diagram to derive a best-fitting absolute magnitude 
M = —19.31 + 5 log 10 (/i/0.7), which was subtracted from the corrected magnitudes 
to give the listed distance moduli. The latter thus depend on the choice of Hubble's 
constant; Astier et al (2006) adopt h = 0.7 as their fiducial value. Although knowledge 
of these parameters is not essential to the determination of cosmological parameters, 
which depend only upon magnitude differences, they are nevertheless useful adjuncts. 
Riess et al (2004) follow a similar procedure, but do not mention values of H and M 
explicitly (see the discussion in their Appendix A); there is no reason to suppose that 
the implict values are the same as those used by Astier et al (2006). Having reduced 
the gold values by 0.14 magnitudes before combining the two lists, it is a matter of taste 
whether we use distance moduli or corrected magnitudes as working variables; we prefer 
the latter, and have reconstructed a table of same by adding —19.31 to each \i B in our 
composite list. In what follows we marginalize over the absolute magnitude M, rather 
than treat it as a fixed parameter. 

3. Cosmological parameters 

With respect to the milliarcsecond radio-sources, our procedure follows that given in 
Jackson (2004). Of the 613 sources mentioned in Section 2, 468 have z > 0.5, with 
0.501 < z < 3.787. These are placed in appropriate redshift bins, and the mean redshift 
and a suitably defined characteristic angular size 9 relating to each bin comprise our 
data points; our working numbers are logs of 9. As in Jackson (2004), we find that 
the best compromise (between many bins containing few objects and a small number 
containing many) is to have 6 bins, which conveniently is a divisor of 468 to give 
78 objects in each. This number of bins is comfortably larger than the number of 
parameters to be determined, and the correspondingly large number of objects within 
in each accurately determines the characteristic angular size. To derive confidence 
regions, we give each point equal weight, and define a corresponding standard deviation: 
a 2 = residual sum-of-squares/ (n—p), relative to the best-fitting curve defined by f2 m , 
and d, where n = 6 is the number of points and p = 3 is the number of fitted parameters; 
here d a characteristic linear size associated with the source population. The definition 
of 'characteristic' in this context is discussed at length in Jackson (2004), where it is 
shown that there is a clear inverse correlation between linear size and absolute radio 
luminosity. In a flux-limited sample those sources with large redshift are intrinsically the 
most powerful, which means that the corresponding size distribution is biased towards 
intrinsically smaller objects; this introduces a weak bias towards open universes, if 
characteristic angular size within each bin is taken to be the simple mean. As a 
countermeasure the characteristic angular size is defined as a median biased towards 
the lower end of the size distribution, namely the boundary between the bottom third 



Legacy data 



7 



and the top two thirds of the latter. With 78 objects in each bin the boundary is defined 
by the 26th point; in fact we take the mean of points 16 to 36 (i.e. 26±10) as a somewhat 
smoother measure than the 26th point alone, but this range is not critical. With this 
definition we find best-fitting values Q m = 0.22, Q A — 1-06 and d = 7.Toh~ l kpc, and 
a standard deviation a = 0.0074 in log^, i.e. 1.7% in 6. To show that these figures 
are robust with respect to binning, we have repeated the calculation for 15 bin sizes, 
from 13 bins of 36 objects to 6 bins of 78 objects, in steps of 3 objects, retaining the 
maximum redshift of 3.787 in each case and using the largest number of bins compatible 
with having no objects with z < 0.5. The outcome is f2 m = 0.24±0.04, = 1.01±0.07, 
and curvature parameter K = 0.25 ± 0.04 (68% confidence limits). These figures are 
consistent with those derived from a confidence region, discussed below. The question of 
robustness with respect to the definition of characteristic angular size will be considered 
at the end of this section. 

Rather than rely too heavily upon a particular binning, we have produced points 
corresponding to bin sizes of 76, 77 and 78 objects, for which the respective standard 
deviations in log 6 are a = 0.0125, 0.0031 and 0.0074; our working points are a composite 
of these three cases. We do not add smaller bins, because this generally means discarding 
valuable data close to z = 0.5. The resulting data points are shown in Table ^ the 
standard deviation is now a = 0.00603. The red plots in Figure H3 shows 68% and 95% 
confidence regions in the £7 m -f2A plane, marginalized over the nuisance parameter d; 
1-dimensional projections give Q m = 0.25 + 0.04/ - 0.03, tt A = 0.97 + 0.09/ - 0.13 and 
K = 0.22 + 0.07/ - 0.10 (68% confidence limits). 

Table 1. Data points for the angular-size/redshift relationship; 9 is in milliarcseconds. 



z 9 


0.6153 


1.4624 


0.8580 


1.2801 


1.1527 


1.1599 


1.4200 


1.1448 


1.8288 


1.1760 


2.5923 


1.2374 



With respect to the acoustic horizon, the key observational point is the position 
of the first peak in the CMB angular spectrum, now located precisely at multipole 
moment / pea k = 220.1 ± 0.8 (Hinshaw et al 2003). We must first specify values of H 
and Q^h 2 , and will adopt the Wilkinson Anisotropy Microwave Probe (WMAP) values 
H = 72 ±5 km sec" 1 Mpc' 1 and Vt h h 2 = 0.024 ±0.001 (68% confidence limits) (Spergel 
et al 2003). Strictly speaking we should not use values which derive from the angular 
spectrum, but these figures are well supported by independent observations, Cepheid 
variables in the case of Ho (Freedman et al 2001; Altavilla et al 2004; Riess et al 
2005), and nucleosynthesis considerations for Q^h 2 (D'Oderico et al 2001; Kirkman et 
al 2003; Pettini and Bowen 2001). If equation (P) is used to calculate the size of the 
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acoustic horizon the biggest source of uncertainty arises from the uncertainty in Hubble's 
constant; ZacH changes by ±2.9% as H changes over the above error range, whereas the 
corresponding figure for Q^h 2 is ±0.28%. Thus our procedure is to fix Q^h 2 at 0.024, set 
^AcH = gU(X) x vr/220, and regard equation (JTJ) as giving h as a function of Q m and ^a- 
Equation (JTJ is solved to give h at each point in the f2 m -f2A plane, where dx(z r ) is the 
corresponding angular-diameter distance; we then define values \ 2 — [{h — 0.72)/0.05] 2 , 
which with 1 degree of freedom gives the appropriate confidence region. Note that the 
error in / pea k is quite negligible in this context. The continuous blue lines in Figure 
El define 95% confidence limits; the fit is exact along the dashed blue line, which has 
slope —1.255 and intersects the flat (magenta dotted) line at Q m = 0.271. The slight 
curvature in the contours at the low-density end of the diagram is because radiation 
still makes a significant contribution to the total density at z r : (pR/p m )z r = 0.45 when 
tt m = 0.1. 

The green plots in Figure HU show 68% and 95% confidence regions for 157±71 = 228 
SNe la, marginalized over absolute magnitude M. In weighting each point we have 
followed exactly the scheme given in Astier et al (2006); each magnitude is assigned 
a variance a 2 + a 2 nt , where a p is the listed photometric error, and ai nt is a fixed error 
attributed to intrinsic variations, which is chosen to give a minimum value of \ 2 equal 
to the expected figure 225. We find a int = 0.105 for our composite sample; we have 
checked the corresponding figure for the full (nearby plus distant) SNLS sample, and 
find ai n t — 0.131, which coincides with the value given in Astier et al (2006). 

A summary of best-fitting cosmological parameters is given in Tables El and El for 
the three basic samples separately and in various combinations. 



Table 2. Best-fitting cosmological parameters, 95% confidence limits; 

MAS=milliarcsecond radio-sources, AcH=acoustic horizon, SN=Type la supernovae. 



Sample 




f2 m 






n A 






K 




MAS 


0.25 + 0.09/- 


-0.06 


0.97- 


-0.15/- 


-0.34 


0.22 + 0.11/- 


0.26 


MAS+AcH 


0.31 - 


-0.08/- 


-0.06 


0.69- 


-0.09/- 


-0.11 


-0.01 


+ 0.04/- 


-0.05 


SN 


0.39- 


-0.18/- 


-0.21 


0.86- 


-0.31/- 


-0.38 


0.25 4 


-0.47/- 


0.57 


SN+AcH 


0.30- 


-0.07/- 


-0.06 


0.69- 


-0.08/- 


-0.09 


-0.01 


+ 0.04/- 


-0.04 


MAS+SN 


0.31 - 


-0.05/- 


-0.04 


0.76- 


-0.12/- 


-0.15 


0.07 4 


-0.12/- 


0.14 


MAS+SN+AcH 


0.30- 


1-0.05/- 


-0.04 


0.69- 


-0.06/- 


-0.07 


0.00 4 


-0.04/- 


0.04 



Figure |U is an alternative presentation of the same information, in which the 
coordinates have been transformed to the deceleration parameter q = Q m /2 — and 
curvature parameter K = Q m — Q\ — 1. The best figures for q are —0.55 ± 0.09 (95%) 
and ±0.05 (68%). Figure 0] is clearly compatible with the conclusion that we are living 
in a spatially flat Universe with accelerating expansion. 

The outstanding question relates to robustness of the MAS figures with respect 
to definition of characteristic angular size. If we use simple means instead of biased 
medians, 68% figures for MAS alone become Q m = 0.28 + 0.06/ - 0.04, tt A = 



Legacy data 



9 



Table 3. Best-fitting cosmological parameters, 68% confidence limits; 
MAS=milliarcsecond radio-sources, AcH=acoustic horizon, SN=Type la supernovae. 

Sample fl m Oa K 

MAS 0.249 + 0.041/- 0.031 0.969 + 0.089/- 0.135 0.218 + 0.066/- 0.103 

MAS+AcH 0.308 + 0.038/- 0.033 0.687 + 0.047/- 0.054 -0.005 + 0.022/- 0.024 

SN 0.389 + 0.096/- 0.104 0.856 + 0.165/- 0.183 0.245 + 0.249/- 0.277 

SN+AcH 0.300 + 0.033/- 0.032 0.694 + 0.040/- 0.042 -0.006 + 0.020/- 0.020 

MAS+SN 0.309 + 0.024/- 0.024 0.763 + 0.065/- 0.071 0.072 + 0.063/- 0.070 

MAS+SN+AcH 0.304 + 0.024/- 0.023 0.693 + 0.034/- 0.035 -0.003 + 0.021/- 0.019 



0.74+0.18/-0.30 and K = 0.02+0.14/-0.25, which represent a significant shift towards 
open universes, but not to an extent which alters the overall qualitative picture. This 
shift illustrates the point made in the first paragraph of this section, and is precisely 
the effect found in Jackson (2004). With respect to composite data sets the shift is of 
virtually no consequence, figures for MAS+SN+AcH becoming Q m = 0.29 + 0.03/ — 0.02, 
tt A = 0.70 + 0.03/ - 0.04 and K = -0.01 + 0.02/ - 0.02. 

4. Evolution 

We do not attempt here to map out the evolution of vacuum energy with cosmic time 
t. The concensus is that large volumes of high-quality data will be needed to do this 
in an unambiguous way (Aldering 2005; Mosoni et al 2006), an order of magnitude 
greater than those available currently. The vacuum equation-of-state is characterised 
by the ratio of vacuum pressure to vacuum energy density, denoted by w(z). Simple 
functional forms have been used to model evolution, for example w = wo + w±z (Cooray 
and Huterer 1999) and w = w + Wiz/{l + z) (Linder 2003); however, Maor et al (2002) 
have shown by means of simulated data and counter-example that such parametrizations 
can produce very misleading results. If w is assumed to be constant then the best-fitting 
value does not necessarily coincide with the time-averaged value. In general there is a 
propensity to produce apparently well-fitting functions w(z) which are distinctly more 
negative than they should be, in some cases significantly less than —1; the latter are 
generally regarded as unphysical (see also Maor et al 2001; Linder 2004; Jonsson et al 
2004; Bassett et al 2004, Jassal et al 2006). There are also non-parametric approaches 
(Chiba and Nakamura 2000; Daly and Djorgovski 2003, 2004; Gerke and Efstathiou 
2002; Huterer and Turner 1999, 2001; Huterer and Starkman 2003; Wang and Freese 
2006) 

As an alternative to the simple functional forms w(z) mentioned above, we shall use 
a proper physical model of dynamical vacuum energy, based upon the simple potential 

V(4) = (8) 

where ujq is the Compton frequency of the corresponding boson. If it turns out that the 
data favour a value ojq > 0, we shall take this as evidence that the Universe is governed 
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by a scalar vacuum energy which is diminishing with time in a reasonably smooth way, 
rather than a cosmological constant, and not as support for the particular potential (JHJ). 
In other words we are using the latter to test the hypothesis that the rate-of-change of 
vacuum energy is zero. In what follows we shall assume that the Universe is spatially 
flat, and that the true cosmological constant is strictly zero. Corresponding to the said 
potential, we have scalar density 

P0 = (0 2 + ^0 2 )/2 (9) 

and pressure 

^ = (0 2 -c^</> 2 )/2. (10) 
The relevant Friedmann equation for the scale factor R[t) is 

4:71 4:71 

R = - yG(p +3^+p m +2p R ) J R = -— G(20 2 -^ 2 2 +p m +2p R ) J R (11) 

and the scalar field is governed by the equation 

4> + 3H<p + u? c (j) = (12) 

where H = R/R is the instantaneous value of Hubble's 'constant'. 

To solve the system of equations (fTTj) to (fT2|). it is convenient to introduce 
dimensionless variables x = R/Rq, y = x = R/ Rq, and 



4nG\ 1 / 2 ± (4kG\ 1 ' 2 

— — 4> V = U = — 



u = [—) v = u = (13) 

where the dot is now the derivative with respect H t, in other words our unit of time is 
the Hubble time. Equations (fTTj) to (fT2"j) are then re-cast in the form of an autonomous 
non-linear dynamical system: 

ii = v v = —3yv/x — UrU 

(14) 

x = y y = -{2v 2 -u 2 c u 2 )-\njx 2 -n K /x\ 

Here ujq is ujc/Hq in the old units (rather than introduce a new symbol we trust that 
context will dictate the appropriate form) and = 8.023 x 10 -5 , corresponding to 
h = 0.72. In terms of the new variables the angular-diameter distance is 

d A ( Z = (I + Z)- 1 — I j-r ■ 15 
H Jt(z) X{t) 

By definition initial conditions (at t = t ) are x(to) = 1 and y(to) = 1; those 
relating to u and v are in part determined in part by a specified value of the scalar 
density parameter fl^ = 87rGp^(to) /3Hq, which using equations © and (|T3*j) becomes 

n <t> = LU 2 c u(t ) 2 + v(t ) 2 . (16) 

However, equation (JT2~j) is second order, and we must specify initial values u(t ) and 
v(to) = u(t ) to give a unique solution, or equivalently current values of the scalar 
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density and pressure. We need to find an extra condition which fixes u(t ) and v(t ) 
unambiguously. To this end, we note that when z > 3 we expect matter or radiation 
to be dominant, and R(t) oc t n , with n = 2/3 or 1/2 respectively; equation (JT2J) then 
becomes 

+ ^0 + ^0 = 0. (17) 

Equation (fTTj) has a regular singular point at t = 0, and has a two-parameter series 
solution, for example 

<f> = A[l- (uct) 2 /6 + ...} + (uj c t) 2 /2 + ...] (18) 



when n = 2/3, and 



B 



= A[l - (uct) 2 /5 + ...] + ^ [1 - (wct)73 + •••] (19) 

when n = 1/2. The problem is that as we integrate backwards (t < to) the last term in 
equations (fTSjl or (fT§j) grows rapidly, and the corresponding vacuum energy can become 
dynamically dominant at critical cosmological epochs, particularly nucleosynthesis, 
recombination and galaxy formation. This must not be allowed to happen, and the 
singular component must be eliminated (Ratra and Peebles 1988). The problem is not 
peculiar to the particular potential (jHJ), but is encountered whenever we use equation 
(|12|) with an explicit potential as a model of vacuum energy; in earlier work the problem 
was avoided by setting initial conditions at t — 0: 0(0) = and 0(0) fine-tuned to 
give the desired current conditions, as in Jackson (1998), Jackson and Dodgson (1998) 
and Weller and Albrecht (2002). There is a possible degeneracy here, in that strictly 
speaking there is an allowable range of non-zero values of |0(O)| << c<j c |0(O)|; however, 
the constraints imposed by the above-mentioned astrophysical are constraints are severe, 
and within the allowable range the the results to be presented here are not noticeably 
dependent on the value of |0(O)|. In the present context we find the above-mentioned 
approach somewhat cumbersome, and have adopted an alternative means of achieving 
the same ends, as follows. 

We are almost certain that the expansion is accelerating, which in this context 
means that is close to its slow-roll value and is falling slowly towards the minimum in 
V(0) at = 0; during this phase we expect to be given by 

*=-g (20) 
which as an initial condition on v would be 

v(t ) = -u 2 c u(t )/3. (21) 
Equations (fTHjl and (|2~T|) together give 

/ O \ 1/2 

«(*>) = 2 / 4/Q • (22) 
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We have tried expressions (}2T]) and as starting values, and find them to be 
wanting. Instead we must retain a degree of freedom, and use 



The procedure is to integrate back to a suitably high fixed redshift (depending on the 
application, for example Z{ = z T = 1088 to determine dA(z T )), and to find Q m (zf) as a 
function of the variable 7. The correct value of 7 is the one which maximizes f2 m (7, zi). 
As an indication of the necessity of this undertaking, we quote figures for a model with 
current values fl m = 0.3, fl^ = 0.7; with u c = we find f2 m (7, 10) = 0.9956, the balance 
comprising scalar energy and radiation; with ujq = 1 we find f2 m (3, 10) = 0.1140, whereas 
the optimization procedure gives fi m (3.59, 10) = 0.9952. In other words the slow-roll 
value 7 = 3 gives a universe in which the scalar component becomes increasingly 
dominant as we move to earlier epochs, whereas for the correct value 7 = 3.59 the 
opposite is the case. 

Figure shows 68% and 95% confidence regions in the uq-VL^ plane, to be discussed 
below. 

5. Conclusions 

In isolation, none of the techniques currently used to constrain cosmological parameters 
give precise numerical answers; at best we get answers to profound qualitative questions, 
relating particularly to acceleration and deceleration. It is only in combination 
that such techniques give meaningful numerical results. Nevertheless, we believe 
that our figures relating to milliarcsecond radio-sources are the most accurate single- 
technique ones to date, fully justifying the expectations expressed in the last paragraph 
of the Introduction; those relating to the joint samples MAS+Ach, SN+AcH and 
MAS+SN+Ach (TablesElandEl) are close to the best currently available (see for example 
Tonry et al 2003; Riess et al 2004; Astier et al 2006). Figures El shows that there 
is no evidence here for dynamical vacuum energy, in that the data prefer the value 
uc = 0. However, the data are compatible with a current rate-of-change corresponding 
to u>c = 1-34 times the Hubble rate, which is thus faster than the latter and must be 
regarded as rapid. If the vacuum energy really is constant, the best we will ever do is 
to place limits on its variation, which is a potential dilemma for physics, because the 
smallest variation would suggest something like a scalar model, whereas true constancy 
would imply something more fundamental. 

Figures |H1 and [7| are graphical summaries of the data used in this investigation. The 
dashed blue curve is for a model with fl m = 0.304, Q\ = 0.693. Mock radio sources 
were generated from supernovae using the angular-diameter distances defined by each 




(23) 



and 




(24) 
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apparent magnitude and a fixed absolute magnitude M = —19.28 + 5 log 10 (/t/0.7), 
and a fixed linear size d = 6.839/i~ 1 pc; M and d are the best-fitting values for the 
samples SN and MAS respectively, when Q m and are fixed at the above values. 
Mock supernovae were generated from radio sources in similar fashion. The mock radio 
source corresponding to the acoustic horizon was generated using d, and an angular 
diameter distance defined by its observed angular size and the value of of /a c h given 
by equations © to @, with h = 0.72, tt h h 2 = 0.024 and fi m = 0.304. With these 
definitions the data points reflect the observational spread. 

Postscript. As this paper was about to be submitted the WMAP three year results 
became available (Spergel et al 2006; Hinshaw et al 2006). We have had a quick look 
at the new data, to assess their impact on the figures presented above. According 
to the scheme adopted above we should use the new best-fitting flat ACDM values 
h = 0.73 ± 0.03, Q^h 2 = 0.0223 and / pea k = 220.7. The main effect of these changes is 
to move the AcH region in Figure El slightly to the left, so that the intersection with 
the flat line, previously at Q m = 0.271, is now at Q m = 0.240. This is in line with the 
WMAP three year flat ACDM values, alone and in combination with a number of other 
astronomical data sets (Spergel et al 2006); the exceptions are WMAP + supernovae (the 
gold subset of Riess et al (2004)) and WMAP+weak lensing data, which continue to 
favour somewhat higher values. Our new figures are summarized in Table HJ with the 
exception of the curvature parameter K, these are virtually the same as those in Table 
El the new values of K indicate a marginally open Universe; the corresponding error 
bars are smaller, because the blue region in Figure El is now narrower, but the conclusion 
K < is not yet statistically significant. 



Table 4. Best-fitting cosmological parameters, 68% confidence limits, taking 
into account the WMAP three year results; MAS=milliarcsecond radio-sources, 
AcH=acoustic horizon, SN=Type la supernovae. 



Sample 




f2 m 












K 




MAS+AcH 
SN+AcH 
MAS+SN+AcH 


0.311- 
0.297- 
0.303- 


-0.037/- 
-0.033/- 
-0.025/- 


0.034 
0.031 
0.022 


0.673- 
0.690- 
0.684- 


1- 0.046/- 
h 0.039/- 
h 0.030/- 


0.050 
0.042 
0.033 


-0.016- 
-0.013- 
-0.013- 


-0.017/- 
-0.012/- 
-0.012/- 


0.017 
0.013 
0.012 
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Old S2290 



Figure 1. Legacy data. Plot of flux density S2290 derived from the Parkes catalogue 
PKS90 (New S2290) against the corresponding figure from Preston et al (1985) (Old 
S2290), for those sources (346) which have measures in both. Values of S2290 derived 
from PKS90 are used in suitable cases when a direct measure is not given in Preston 
et al (1985). 
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Figure 2. Multipole moment /p ea k corresponding to the first acoustic peak in the 
CMB angular spectrum; blue line is the analytical curve according to equations (J), 
© and ©, with z T = 1088, h = 0.7 and n h = 0.04; data points are from CMBFAST; 
all spatially flat: £l\ = 1 — O m . The yellow horizontal line is the WMAP figure 

^peak = 220. 




Figure 3. Confidence regions in the f2 m -f2A plane, 95% and 68%; red=milliarcsecond 
radio-sources, blue=acoustic horizon, grccn=Typc la supernovae; the filled areas are 
joint confidence regions for all three; the acoustic horizon fit is exact along the dashed 
blue line. 
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Figure 4. Confidence regions in the qo~K plane, 95% and 68%; red=milliarcsecond 
radio-sources, blue=acoustic horizon, green=Type la supernovae; the filled areas are 
joint confidence regions for all three; the dashed blue lines are now 68% confidence 
limits. 
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Figure 5. Confidence regions in the ojq-H,^ plane, 95% and 68%, where ujq 
is Compton frequency in Hubble units; as determined by the three data sets: 
red=milliarcsecond radio-sources, blue=acoustic horizon, green=Type la supernovae; 
the filled areas are joint confidence regions for all three. Flat universes with evolving 



scalar vacuum energy governed by a potential V(<j>) 



2 /2. 
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Figure 6. Angular-diameter/redshift diagram; each lower blue star is a bin containing 
c. 77 milliarcsecond radio-sources (see text); the error bars are too small to be seen on 
this scale, amounting to ±0.006 in log (9. Each cyan point corresponds to a supernova 
transformed into a mock milliarcsecond source. The top blue star represents the 
acoustic horizon; the error bar is ±0.0120 in log 6, corresponding to h = 0.72 ± 0.05. 
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Figure 7. Magnitude/redshift diagram; each cyan point is a supernova. Each blue 
star represents a bin containing c. 77 milliarcsecond radio-sources (see text); the error 
bars are too small to be seen on this scale, amounting to ±0.03 magnitudes. 
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